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We study bond percolation of N non-interacting Gaussian polymers of £ segments on a 2D 
square lattice of size L with reflecting boundaries. Through simulations, we find the fraction of 
configurations displaying no connected cluster which span from one edge to the opposite edge. 
From this fraction, we define a critical segment density Pc(£) and the associated critical fraction of 
occupied bonds Pc(£), so that they can be identified as the percolation threshold in the L — » oo 
limit. Whereas Pc(£) is found to decrease monotonically with £ for a wide range of polymer lengths, 
Pc{£) is non-monotonic. We give physical arguments for this intriguing behavior in terms of the 
competing effects of multiple bond occupancies and polymerization. 
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Studied for many decades, the percolation problem is 
rich and venerable [1] . To this day, interesting new issues 
continue to surface, for example, in the context of perco- 
lation of circles [2], ellipses [3], random surfaces [4] and 
polymers [5] . Motivated by recent experiments on trans- 
port of small gas molecules across a polymer membrane 
[6], we investigate a novel aspect, namely, bond perco- 
lation with a specific correlation between the bonds. In 
particular, we consider a standard model for polymers: a 
chain of I segments linked to form a random walk (RW) 
in space. If, instead of continuous space, we use discrete 
lattices, it is natural to place the segments on the bonds. 
With certain considerations for our physical system in 
mind, we study ordinary (i.e., non self-avoiding) RW's. 
Clearly, the locations of the bonds occupied by one poly- 
mer are strongly correlated, except for the I = 1 case. If 
we place N such polymers randomly on a finite lattice 
and ask when the occupied bonds percolate, we face a 
"correlated-bond" percolation problem. 

Before defining our problem in detail, let us devote a 
paragraph to a typical experiment. First, a thin poly- 
mer membrane is formed by quenching a polymer melt 
to a temperature far below its glass transition. To study 
permeation by small gas molecules, a pressure gradient 
is set up across the membrane and the gas current, in an 
effective steady state, is measured. The focus of the ex- 
periments in [6] , namely, how aging the membrane affects 
such transport properties, serves only as a motivation for 
our model and investigation here. Detailed studies of gas 
transport will be published elsewhere. 

For simplicity, we introduce a two-dimensional (2D) 
version of this system. As the polymers originated from 
a rubbery state, we feel justified to neglect the complica- 
tions due to self-avoidance; further, the polymer matrix 
is essentially static on the time scales associated with 
gas transport. Therefore, we represent each polymer as 
a simple RW, on a square lattice of size L 2 (cf. Fig. 1). 
To model a finite system with real boundaries, we im- 



pose brick wall boundary conditions on these walks. To 
be specific, we start a "walker" at a randomly chosen site, 
with each step to a nearest-neighbor site (a bond) mod- 
eling a segment of the polymer. When a walker arrives 
at a boundary site, the probability for reversing direction 
is 1/2 (instead of 1/4 in the bulk). In this way, a bond 
represents a "segment" of the polymer of the persistence 
length, rather than a physical monomer. Note that, since 
the RW's are not self-avoiding, each bond may be mul- 
tiply occupied and we will use the term "m-occupied" 
to signify a bond carrying m segments. Meanwhile, the 
gas molecules are located on the cells of our lattice and 
diffusion is modeled by hopping to a nearest-neighbor 
cell, across a bond. The rate for a hop across an m- 
occupied bond will be controlled by an activation barrier: 
exp (— me/ksT), where T is the temperature and e is an 
energy associated with jumping over a single segment. 
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FIG. 1. Gas molecules diffusing from left to right across a 
10 x 10 lattice, carrying N = 3 polymers of length I = 20. 

As the polymers are static in this model, we are deal- 
ing with driven diffusion in a quenched random medium. 
Transport through static disordered media with random, 
energy barriers has been extensively studied in the past 
few decades. Unlike many previous studies [7] of this 
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type, the energy barriers in our medium possess strong 
spatial correlations since the segments of a single RW 
are highly correlated. In this letter, we focus on a much 
simpler question: In the low-T, high-e limit, what is the 
probability that there will be any current through the 
membrane? In this limit, crossing any occupied bond is 
forbidden. Thus, for a 2D system, the quantity of in- 
terest is precisely the probability for the occupied bonds 
to form a connected cluster that spans the system in the 
direction transverse to the pressure gradient, i.e., the per- 
colation probability in one of the directions. Of course, 
the picture is more complex for a 3D system, on which 
we will comment at the end. 
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FIG. 2. The probability distribution V{p;£,N,L) for the 
fraction of occupied bonds p for p = 1, I = 256, and different 
system sizes. Inset: The mean value p vs. polymer length t 
for p = 1/2 and p = 1. The RMS deviation for all data points 
is of the order of 10 -3 . 

Let us briefly comment on the differences between our 
focus and earlier studies which addressed this problem 
[5,8]. Previously, all polymers were modeled as self- 
avoiding walks (SAW). Clearly, these pose more chal- 
lenging combinatorial questions than our RW's, so that 
analytic progress is essentially impossible. In another re- 
spect, however, SAW's are simpler: Since bonds are at 
most singly occupied, a given density of polymers maps 
directly into a delta distribution for the probability of 
bond occupation. Further, only relatively short poly- 
mers (£ < 16 in [5]) were considered, focusing on singular 
properties near the percolation threshold. Thus, finite 
size scaling techniques were invoked to draw conclusions 
in the thermodynamic limit. By contrast, the samples in 
permeation experiments are quite thin so that £ and L 
can become comparable. As a result, we will need to pay 
attention to finite system sizes here. To the best of our 
knowledge, ours is the first attempt to study percolation 
of polymers on a lattice with multiple bond occupancies, 
in this particular scaling limit. 

Turning to the specifics of our problem, we place the 
pressure gradient in the x-direction, so that the presence 



of a particle current translates into the absence of a con- 
nected cluster of occupied bonds spanning the boundaries 
y = 1,L. Our question is now well posed: Of all systems 
with N random walks (polymers) of £ steps (segments) , 
placed randomly on a L 2 lattice, what is /, the fraction 
which contains no cluster of connected segments span- 
ning the system in one of the directions? 

Clearly, in a gas permeation experiment, / is just the 
fraction of 'open' configurations, allowing nonzero cur- 
rent. Let us emphasize that, due to multiple occupancies, 
there is no simple, 1-1 relation between p, the fraction of 
occupied bonds, and the segment density, p = N£/(2L 2 ) 
(proportional to the mass density of the sample). In con- 
trast, for the SAW's in [5], p is exactly p. Here, we 
have a distribution of p's, a quantity to be denoted by 
V(p;£, N, L). Even in the I = 1 case - for which single 
segments are randomly placed and we can easily find the 
probability for a bond to be unoccupied: [1 — 1/(2L 2 )] N 
- the entire distribution is nontrivial [9]. For arbitrary £, 
we are not aware of an analytic expression for V . How- 
ever, for large N and L, with fixed £, we may expect V to 
be sharply peaked, as illustrated in Fig. 2. In all cases, 
we may define the average p = ^2 pP and use it as a con- 
trol parameter instead of p. Of course, their relationship 
is of interest. In general, we have no analytic expression 
for p (p; £, L) though we are able to find good estimates in 
some regimes [10]. To demonstrate this nontrivial depen- 
dence, we show in the inset of Fig. 2 simulation results 
for p [p; £, L) vs £ at two different densities and L = 128. 
We find that increasing I at fixed p decreases the fraction 
of occupied bonds, since longer polymers tend to overlap 
more. Note that neither the details of V nor its average 
p are needed for determining the main quantity of inter- 
est: / (£, N, L). We raise these issues related to p purely 
for those familiar with the standard question of ordinary 
bond percolation: If a fraction p of the bonds is occupied, 
what is the probability for the system to percolate? 




FIG. 3. The percolation probability f(£, N, L) vs. the seg- 
ment density p for fixed I = 256 and several L. 
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To study /, we perform simulations for various p G 
[0.25, 2.0] and L G [32, 256]. With a full range of £'s and 
iV's, we find not only the typical percolation transition, 
but also seemingly contradictory effects due to the corre- 
lations. Focusing on the first aspect, we consider the data 
for fixed £ but increasing L. For example, in Fig. 3, we 
plot / against p, with £ = 256 and L = 32, 64, 128, 256. 
Not surprisingly, / is essentially unity for small p and zero 
for large p. In the L — > oo limit, we expect / to become a 
step function, vanishing for p > p c , i.e., the "percolation 
threshold." For our finite systems, we see the approach 
to this behavior, namely, f(p) falling off more steeply 
as L increases. As an estimate for p Cl we construct the 
(average) local gradients, — <9 p /(p), which we expect to 
approach 5 (p — p c ) with increasing L. Exploiting ~d p f 
as a "distribution" over p G [0, oo], we compute the av- 
erage — J °° pdpf and define it as the "critical" p^ (£) for 
each L. For the above L's, these are respectively, 0.849, 
0.754, 0.711, and 0.688. A naive extrapolation (linear in 
1/L) leads to pf{£ = 256) ~ 0.67. To express this result 
in the usual language of percolation, we find p (£, p; L) 
from the numerical data and define the "critical" p via 
PcW = v{^iP = PcWi^)' Extrapolating to L — > oo as 
before yields the percolation threshold p/?{l = 256) ~ 
. 295 (data for L = 32, 64 and 128 shown in Fig. 4). This 
value is consistent with the £ — > oo projected value for 
the SAW's [5]: 0.349±0.001 (based on extrapolation from 
£ < 16). A more quantitative comparison must wait for 
better statistics and more sophisticated analyses (e.g., 
finite size scaling) and will be reported at a later time. 
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FIG. 4. The effective percolation threshold p% (I) vs. 1/L 
for various polymer lengths I. 

Next, we turn to a second, very natural question: If 
we fix the density of segments in a system, how will poly- 
merizing them, into chains of increasing length £, affect 
f(£, N, L)l Here, we discover a rather intriguing aspect 
of this problem, illustrated in Fig. 5: For small densi- 
ties, f(£, N, L) decreases with £, but displays the oppo- 
site trend for high densities - increasing with £. For in- 
sight into these seemingly contradictory findings, we pro- 
vide some intuitive arguments. At low densities, linking 



isolated segments together into longer chains clearly en- 
hances the chances of spanning. As an extreme example, 
consider L segments so that only one configuration can 
span the lattice. For £ = 1, at the least, all of them must 
have the same orientation and lie in the same column, 
so that the probability of spanning is o((2L)~ L ). Yet, 
if they are linked to form a single polymer, this prob- 
ability rises to 0(4~ L ). In contrast, at high densities 
correlations have the opposite effect: if randomly dis- 
tributed segments percolate, linking them together into 
chains tends to localize them and enhances overlaps. As 
a result, percolating paths can be broken, and / (£, N, L) 
increases. 
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FIG. 5. The percolation probability f(£, N, L) for fixed 
density p — N£/(2L 2 ) is plotted against £ for system size 
L = 128. 



0.76 
0.75 
0.74 



g 0.73 
^ 0.72 
0.71 
0.70 






1 1 

L = 64 + 






1 = 128 -X- 






\\ / 












V'"' 






1 1 





1 



10 



I 



FIG. 6. The critical segment density p^(£) vs. £ at the per- 
colation threshold for two system sizes. For L — » oo, this plot 
may be viewed as a phase diagram: above the line, the poly- 
mers form a percolating path, while none exists below. The 
data points correspond to £ — 4™ where n = (0, 1, 2, 3, 4). 

These competing effects of polymerization induce a 
dramatic ^-dependence in the critical density p^. As 
illustrated in Fig. 6, there appear to be three distinct 
regimes, with two extremal points for Pc(£). Let us pro- 
vide some intuitive arguments for these properties. Fo- 
cusing on the first regime (£ < 4), consider £ = 1, where 
our data suggest a critical density of p^°(l) — 0.69 (in ex- 
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cellent agreement with the theoretical In 2 ~ 0.693). If we 
now generate the same segment density (0.69) with poly- 
mers of length 2, the second segment has a probability of 
1/4 to overlap with the first one, so that this density will 
not suffice for percolation. Thus, p^(2) must be higher 
than (1), in agreement with our observations^ (2) rj 
0.74 in simulations). This effect continues to dominate 
up to a certain value £* ; from the data shown in Fig. 6, we 
may conclude 4 < £* < 16. Our data indicate only a very 
weak L-dependence in this regime, so that we are confi- 
dent that the L — > oo limit is essentially reached. Next, 
we consider the third regime: polymers with £ > L (the 
last two data points in Fig. 6). Here, each polymer al- 
ready occupies a significant fraction of the system. Thus, 
adding another increases p noticeably, but changes bond 
occupations only marginally. In this sense, a dispropor- 
tionately larger number of polymers is required to span 
the system, leading again to increasing p^ with I. Since 
this regime is dominated by the finite size of the system, 
we expect it to shift to larger £ for bigger L's. Sandwiched 
between these two regimes, we observe the opposite ef- 
fect, due to polymerization: After reaching a maximum 
at £*, the critical density decreases. In this regime, for 
reasons yet to be understood, the argument we presented 
above applies again: tieing segments together generates 
"longer" objects that percolate more easily. We conjec- 
ture that this second regime extends to infinity in the 
thermodynamic limit, so that p^°(£) exhibits only a sin- 
gle extremal point at £*, with 4 < £* < 16. 

In contrast, when expressed in terms of the bond oc- 
cupancy, we observe no surprises: p^ (£) is monotonically 
decreasing. Since we also expect a monotonic depen- 
dence of p on p (for fixed £) , it is all the more remarkable 
that the effects of correlations built into the segments by 
polymerization are so subtle to produce non-monotonic 
behavior in Pc{£)- 

In summary, we have investigated f(£, N, L), the prob- 
ability that N Gaussian polymers of £ segments, ran- 
domly distributed on the bonds of a square lattice, do 
not span anLxL system. Since the polymers are sim- 
ple random walks, the bonds are frequently occupied by 
more than one segment, leading to a nontrivial corre- 
spondence between segment density, p = N£/(2L 2 ), and 
bond-occupation probability, p. For a given £, a sharp 
percolation threshold emerges as L — > oo. For finite L, 
we define "critical" values p^ (£) and p^ (£). Measuring 
them by simulations, we study their behavior for various 
£. We find that the subtle interplay of polymerization 
and multiple occupancies leads to dramatic effects, espe- 
cially for Pc{£). Several limits, e.g., the extreme case of 
f(£ = 2pL 2 , 1, L), are analytically accessible and will be 
reported elsewhere [10]. 

It is worthwhile to devote a few lines to comment on 
the experimental implications of our work. Although our 



simulations were restricted to two dimensions, we believe 
that our principal result, i.e., the two opposite effects of 
polymerization on percolation probability (lowering p c 
on one hand, while making more free bonds available 
via multiple bond occupancies on the other) should be 
present in higher dimensions also. Percolation-related 
transitions in permeation properties have been observed 
experimentally, eg. in [11] where a jump in the gas cur- 
rent is seen when the volume fraction of the flexible com- 
ponent of the polymer is increased beyond a threshold 
value. It would be interesting if the non-monotonic vari- 
ation of the critical density with polymer length could be 
tested experimentally. 
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